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We discuss general aspects of the possibility of Bose condensation of diquark pairs 
in systems of dense matter. Lattice field theory simulations are presented for model 
four-fermion theories which are expected to manifest the phenomenon, and results 
from measurements of both diquark two-point and one-point functions presented. 
Whilst inital results are promising, there remain systematic effects needing to be 
understood. 



There has been much recent interest in the possibility of Bose condensates 
f*^ ■ formed from quark pairs close to the Fermi surface in dense matter; such con- 

densates, which carry non-zero baryon number, can. lead to exotic scenarios for 
dynamical breaking of the color gauge symmetry EJ. So far most work in this 
direction has relied on effective field theory descriptions of the strong interac- 
tion. For any problem invloving the QCD ground state it would be nice to have 
a calculation which in principle involves no uncontrolled approximations; this 
is therefore a natural area for the application of the methods of lattice gauge 
O i' theory. In this talk we describe the first steps we have made in this direction; 

our initial results have already appeared 0. 

When thinking about possible diquark condensates (qq) ^ which might 
form there are several issues to consider. Firstly, is the condensate wavefunc- 
tion gauge invariant? In the QCD scenarios the condensate is formed from a 
3 <g> 3 of the color gauge group, which necessarily breaks the local symmetry, 
leading to a dynamical Higgs mechanism making some or all of the gluons 
massive. This is the phenomenon of color superconductivity. In other models, 
relevant to this discussion, another possibility is that the condensate wavefunc- 
tion is gauge invariant; in this case condensate formation breaks only global 
symmetries, and yields a superfJuid state. Next, is the condensate necessarily 
a spacetime scalar? We normally expect condensates to share the spacetime 
quantum numbers of the perturbative vacuum, but in systems of dense matter 
there is a preferred rest frame, making it possible to consider rotationally non- 
invariant condensates El. Another possibility which should not be excluded a 
priori is that the condensate spontaneously breaks parity invariancen. Finally, 
the most crucial aspect of the condensate wavefunction is that it respects the 
Pauli Exclusion Principle, ie. it must be antisymmetric under exchange of all 
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possible quantum numbers between the quarks. It is this requirement which 
makes the ground state_so exquisitely sensitive to the number of light quark 
flavors present in QCD EJ; more generally it may also result in a sensitivity to 
which representation of the gauge group is carried by the quarks. 

The most generic symmetry we expect to be be broken by a diquark con- 
densate is the global U(l)y of baryon number: 

q^e la q ; q^qe~ la . (1) 

One issue we should address is that spontaneous breaking of a global vectorlike 
symmetry like this is usually forbidden by the Vafa-Witten theorem □. It is 
convenient to classify the models in which diquark condensation might occur by 
the various escape clauses offered by the theorem. The theorem does not apply 
if the path integral measure is not positive definite. This is the case for QCD 
with chemical potential fi ^ 0, since for this case the fermion determinant is 
complex. Unfortunately it is precisely this feature which makes Monte Carlo 
simulation extremely difficult. The theorem also fails if the theory includes 
a Yukawa coupling to a scalar degree of freedom. Two model field theories 
satisfying this criterion are the Gross-Neveu (GN) model in 2+1 dimensions, 
and the Nambu - Jona-Lasinio (NJL) model in 3+1 dimensions, and these will 
be the main focus of this talk. The final escape route is if baryon number 
is part of some larger non-vectorlike symmetry. An example is SU(2) lattice 
gauge theory where the global symmetry group enlarges from U(l)y <g> U(l)^ 
to U(2) due to the pseudoreal nature of the 2 representation. Simulations of 
SU(2) lattice gauge theory with [i + 1 will be discussed separately!! 

The GN and NJL models, which are essentially identical apart from the 
number of spatial dimensions, are relativistic generalisations of the model orig- 
inally considered in the BCS mechanism for superconductivity. In continuum 
notation the Lagrangian may be written: 

C = $tf + m)V - 9 2 [(^) 2 - OA75TV0 2 ]- (2) 

The four-fermi terms can be replaced by Yukawa couplings to scalar and pseu- 
doscalar auxiliary fields. The lattice transcription of (||) is discussed in detail 
elsewherea. In addition to the U(l) baryon number symmetry (uh, the model 
has a SU(2) L (g> SU(2) R axial symmetry in the chiral limit m — > 0: 

ip L i-> Uip L , ip L h-» i) L U ] ; ip R h-> Vi> R , 4> R h-> -ipRV\ (3) 

with U, V independent SU(2) matrices. 

The most relevant features of the model for our purposes are that for strong 
coupling the SU(2) L (g) 811(2)^ symmetry spontaneously breaks to SU(2) y by 
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Figure 1: 

formation of a chiral condensate (ipip). The spectrum in the broken phase 
contains both "baryons" , namely the elementary fermions which now have a 
dynamically generated mass, and "mesons" , namely ipip composites, which 
include 3 Goldstone pions. It turns out that for the 2+1 dimensional GN case 
the model has an interacting continuum limit at the critical coupling required 
for symmetry breaking. Crucially, the model can be formulated on a lattice 
and simulated for \i ^ 0B; it is found that in the broken phase a first order 
chiral symmetry restoring transition occurs for some critical /x c Q. 

We have consideredn the possible formation of a diquark condensate with 
wavefunction 



qq = ip tr C'f 5 (g) t 2 <8> r 2 ip; 



(4) 



the operators in the tensor product denote that the diquark is a spacetime 
scalar, and antisymmetric in both implicit (due to lattice fermion doubling) 
and explicit flavor indices (note that C75 is also antisymmetric). Condensation 
of dj) spontaneously breaks baryon number (|l|) but not axial (Q) symmetry. 
We have examined two possible signals for (qq) ^ 0. 

For large spatial separation the diquark pair propagator is proportional to 
the square of the condensate by the cluster property: 

G(x) = {qq(0)qq(x)) = {qq(0)qq(x)) c + (qq)(qq) => lim G(x) - \(qq)\ 2 . (5) 

x — >oo 

A non-zero condensate should therefore reveal itself as a plateau in the large-t 
behaviour of the diquark timeslice propagator. In Fig. 1 we show the behaviour 
of G(t) from simulations of the GN model on a 16 2 x 40 lattice, performed with a 
value fj, = 0.8 > [x c , so that chiral symmetry is restored. There is clear evidence 
for a stable plateau, especially when compared to G(t) for non-interacting 
fermions, also shown with closed symbols. The square root of the plateau 
height is plotted, together with (ipip), as a function of [i in Fig. 2. We observe 
a large increase in the (qq) signal going from the chirally broken phase into the 



symmetric phase, consistent with the notion that the condensates "compete'a. 
We also observe a small but non-zero (qq) in the low density phase, and a 
parity violating pseudoscalar condensate in the high density phase, although 
probably both signals are finite volume artifacts. 

Unfortunately the height of the plateau remains roughly constant as the 
spatial volume of the lattice is increased, whereas naively we expect it to be 
extensive. It is therefore not clear whether a true condensation is occurring, or 
whether there are unexplained finite volume effects. A similar behaviour is seen 
in 3+1 dimensional simulatiojis, suggesting that this is not a low dimensional 
artifact as originally thought a. 

To attempt to clarify these ambiguities we have now performed direct mea- 
surements of {qq} in the simulations. This involves including explicit diquark 
source terms in the action and using a Gor'kov representation, ie: 

S ferm = $M4> + J V tr r 2 V + #r 2 ^ = $, V tr ) ( _ IZtr \ M * ( ' ' ' 



W tr jr 2 J \4> 
= & r A\jj)*; (6) 

Z[jj] = {Pf(A[jj])). (7) 

The diquark condensate is now defined by 
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which is straightforward to implement. Our results are "quenched" in the sense 
that we have included j ' ^ in the measurement routines but not as yet in 
the update algorithm. We show (qq(j)) for various /i for the GN model in Fig. 
3 and for the NJL model in Fig. 4. In the low density chirally broken phase 
the signal extrapolates linearly to zero as j — > 0; in the high density phase the 
signal is larger and considerably less linear, although still vanishing in the zero 
source limit. Preliminary indications are that the curvature of (qq(j)) is still 
more marked in the 3+1 dimensional simulation. We have found evidence for 
small finite volume effects, but the question of whether a diquark condensate 
forms, as revealed by a non-zero intercept in the thermodynamic limit, is still 
open. 

Whilst we have still not found unambiguous evidence for diquark conden- 
sation, it is clear something interesting is happening in the dense phase, as 
revealed both by the long range timelike order in the behaviour of the two- 
point function, and the non-linear behaviour of the one-point function with 
j. Systematic effects due to the influence of Goldstone modes, the relatively 
small number of states close to the Fermi surface on a finite lattice in a small 
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number of dimensions, and the quenched nature of the (qq) measurement may 
all need to be understood. Interestingly, evidence for diquark condensation is 
much more compelling in simulations of SU(2) gauge theory!. It is important, 
however, to understand how a BCS-like mechanism manifests itself in simple 
four-fermion models before QCD can be tackled (recall that in this case the 
condensate is not even gauge invariant). In future work we plan to explore 
the spectroscopy of the model using the Gor'kov representation, and hopefully 
measure the gap, which is after all the quantity closest to physics. 
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